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^\ , We outline two approaches for studying the clectrowcak phase transition in the 

Q\ ■ framework of the four-dimensional SU(2) Higgs model on a lattice. The first one is 

based on a combination of variational estimates for the free energy and a solution 
of the resulting gap equations by means of dynamical linked cluster expansions. 
In the second approach we first indicate the derivation of an effective scalar field 
theory upon integration over the gauge fields. The phase structure of the resulting 
»-^H, scalar theory is then investigated by means of generalized linked cluster expansions. 

We present results for the critical line and the critical endpoint. 
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1 Introduction 

We determine the critical line and critical endpoint of the electroweak tran- 
sition in the framework of the SU(2) Higgs model on the lattice in four di- 
mensions. Usually the electroweak transition in the SU(2) Higgs model is 
studied with Monte Carlo simulations either in the four-dimensional model □, 
cl or -after dimensional reduction- in an effective three-dimensional model of 
the same formB. The time extensions in such Monte Carlo simulations is re- 
stricted to small values such as N T = 2,3 unless one uses asymmetric lattices, 
and A is of O(10 -4 ). In our approach we first want to integrate out the gauge 
field degrees of freedom and then study the phase structure of the resulting 
effective scalar model with Linked Cluster Expansions (LCEs) that have been 
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developed to a powerful tool by extending the expansions to a high order in 
the expansion parameter and to a finite volume Q,U. 

2 Critical Line from Variational Estimates and DLCEs 

In the first part we describe an analytic estimate of the critical line k c {\) 
for a given gauge coupling (3. The estimate is based on a combination of 
variational estimates for the free energy of the SU(2) Higgs model in 4D and 
Dynamical Linked Cluster Expansions (DLCEs) for correlators in the set of 
variational equations. Already in this first estimate the localization of the 
critical line agrees quite well with high precision Monte Carlo results, since 
even this estimate became rather complex after all. 

Variational estimates for the free energy are based on the convexity of the 
cxp-function and the positivity of the measure. Under these conditions we 
obtain an uncquality between partition functions which is of the form 

ex P (-V7) = Z > Z V E-exv<-(S-SvE(0)>VE = exp-V/(C) (1) 

with the following notations. V is the D-dimensional volume, / is the free 
energy density of the system described by the partition function Z, in our case 
Z is the partition function of the SU(2) Higgs model in D=4, Sve (Zve) is the 
auxiliary action (partition function), respectively, that depends on a generic 
set of variational parmeters £, and / is the trial free energy^ density which 
should be optimized as a function of ( in order to minimize / — / > 0. The 
generic set of variational equations is then given as 

dj(0\ c= ^=o d 2 J{()\ c= ~>o (2) 

An equality sign in the second uncquality determines the critical temperature 
T c (coupling k c ) in case of a second order transition and gives a lower (upper) 
bound on T c , (re c ) in case of a first order transition. 

Next we have to find a good choice for Sve in case of the SU(2) Higgs 
model. Since we know that the three-dimensional spatial part of the gauge- 
Higgs interaction is supposed to contain the non-perturbative degrees of free- 
dom that drive the Higgs transition, in our final choice for Sve, we treat the 
physics along three-dimensional spatial hypersurfaces as accurately as possible, 
but implement a factorization along the time direction so that 

Zv E {(link, (cube, = Z c ° be ((cube, () " Z H ° lk 3 ((link) (3) 

in which Zu n f. is an exactly solvable one-link partition function depending on 
a variational parameter (u n k- Lq is the extension in time direction, V3 the 



3-dimensional volume. The action of the partition function Z cu }, e for the 3-D 
hypersurfaces is given as S C ube = S with 

r 3 1 

S = - J^ 2(tTi-^x) + Q ube TYU(x;fi)) + KY,^(&(x)U{x;fiMx + fi)) 

x£A 3 L n=l i 

(4) 
depending on variational parameters Ccube and £. Note that in contrast to 
usual LCEs the hopping term in (0) does depend on U, and U has its own 
dynamics governed by the second term in S cu b e - Thus a Taylor expansion of 
In Z cu b e in powers of k about k = will lead to a more general type of expan- 
sion, since the familiar LCE only apply to frozen [/-dynamics. We will need 
such expansions for derived quantities from In Z cu b e that occur in the set of 
variational equations. The detailed form of the variational equations is given 
intj. Here we only state that they depend on connected n-point functions in 
$s and Us up to n = 6, evaluated w.r.t. < ■ > cu be- It is these expectation 
values that we evaluate as power series in n with DLCEs that have been in- 
troduced and systematized inQ. The number of graphs that contribute up to 
and including 0(k 4 ) is of the order of several hundred. Finally we solve the 
variational equations for Qmk, (cube, S, and k c as a power series in k. 

As results for n c we obtain for example n c = 0.1282(1) (0.1281(1)) for 
A = 5yl0~ 4 (5.M0- 4 ),and/? = 8ascomparedtoK c = 0.12887(1) (0.12852(2)) 
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fromLJ, respectively. The good quantitative agreement of low order DLCEs (in- 
cluding 0(k 4 )) with high precision Monte Carlo results appears less surprising 
in view of the number of graphs that contribute to the series up to that order. 

3 Critical Endpoint of the Electroweak Phase Transition 

Starting from an SU(2)-Higgs model in 4D we have derivedB an effective scalar 
model in the following steps. First we absorb the angular part of the U(2) Higgs 
field in a gauge transformation. The remaining Higgs degrees of freedom are 
then given by scalar fields p(x), p(x) > 0, so that Tr$t$ = p 2 . Integrating 
upon the [/-dependent part of the original action leads to an effective action 
W{p). W(p) is the free energy of the gauge fields in a scalar background field. 
Since we are interested in an effective model that describes the phase structure 
in the vicinity of the critical endpoint with long range correlations, we perform 
a small momentum expansion of gauge field correlators about (lattice) momen- 
tum p = that appear as coefficients in W{p). It is justified if p/m g i ue << 1 
and p/T c << 1 (in a continuum language) with m g i ue denoting the mass gap 
in the pure gauge sector. The small p-expansion of the gauge field correlators 
induces an expansion of the action for the scalar fields with leading ultralocal 



and local terms. Non-local terms are suppressed with higher momenta p. The 
p = term of W(p), W^°\ is calculated in a Monte Caxlo simulation for a 
pure gauge theory in a constant background field r = p(x)u. The Monte Carlo 
results have been compared with results from a small field and large field ex- 
pansion of W in terms of p. It turns out that the final value for the critical 
coupling k c sensitively depends on the shape of W^(nr 2 ) for small values of 
/■sr 2 , but is rather insensitive to the precise form for large arguments. Thus 
we improve on the small field range by including contributions of the order p 2 
to the first term in a small field expansion, but neglect p 2 and higher order 
terms of a large p-expansion, the leading conribution of which is absorbed in 
W°(Kr 2 ). The final form of the effective scalar model S e ff(z) in terms of scalar 
fields z = z(x) = p 2 (x) is then given by 

%(z) ^LA(z-l) 2 +lfW( K z;^ + ^z 2 ^^ z (x)z(^)). (5) 

Although this model looks at most quadratic in z, it is non-Gaussian because of 
the term W^°) resulting from the Monte Carlo determination of the p = 0-part 
of the effective action W and z > 0. For given 0, W^ = W^(kz) is known in 
the form of a table for ~ 500 values of kz € [0, 25]. Note that the model is no 
longer Z(2) symmetric as the original SU(2) Higgs model was in terms of <f>. 
The Z(2) symmetry will only effectively be restored at the critical endpoint. 
The missing Z(2) symmetry is_also responsible for the generalization of LCEs 
to LCEs in an "external field" a. In terms of z the hopping term has its usual 
form, but it should be further noticed that now also two parts of the ultralocal 
part of the action do depend on k. It is possible to account for this additional k- 
dependence by a renormalization of the ultralocal vertices in the LCE graphs: 
the former constant coefficients in the series for the susceptibilities now are 
expanded themselves as power series in k. The phase structure of this model 
is next studied with an effective potential with coefficients calculated to a high 
order in the expansion parameter of the series. 

3.1 A c from an Effective Potential with LCEs 

We approximate the partition function Z(J) of the original SU(2) Higgs model 
in 4D by Z(J) = JVz exp—S e ff(z) + J-z . The effective action T(z) is 
then derived in the standard way by a Legendre transformation of InZ(J) 
with z = dh\Z(J)/dJ . Evaluating T(z) for constant field configurations z 
(corresponding to constant currents J), we obtain the effective potential V e ff. 
Next we express the coefficients of V e ff in terms of quantities that are 
directly available in a high order linked cluster expansion. These are the sus- 



ceptibilities that are 1PI in the LCE- graphical sense and are denoted by Xn PI ■ 
The Xn PI > n — 2, ... ,6, are calculated up to and including 0(k 16 ). If we ex- 
pand the effective potential in fluctuations about the vacuum expectation value 
z := d\nZ(J)/dJ\j = o, the qualitative form looks like 

V e ff(x) = a 2 x 2 + a 3 x 3 + a 4 x i + 0{x 5 ) (6) 

with x = z ~z and coefficients Oj, i = 2, 3, 4 that depend on /?, A, k in a very 
implicit way. The phase structure of the effective scalar model with action S e ff 
is derived by scanning the sign of D := (3a 3 /8a4) 2 — (02/204) as a function of «;. 
A first order transition is indicated by D > for a certain range of ks between 
the occurrence and disappearance of the metastable second minimum in the 
symmetric and broken phase, respectively. In particular the critical endpoint 
(A c , k c ) of the first order transition line for fixed (3 shows up as D = and a 2 = 
= a 3 , or, in terms of the x 1PI s, Ak 2 x\ pi {k, A; j3) = 1 and xl PI { K , A; /3) = 0. 
If we extrapolate the results for A C (L), K C (L), L — 7, .., 16, to infinite order of 
truncation, A C (L) = A c (oo) + const 1/L, we obtain A c (oo) = 0.0032(1) and, 
analogously, k c (oo) = 0.1447(1) . An inclusion of x\ PI an d x@ P7 -terms in the 
effective potential shows that the neglected coefficients in (^) are suppressed 
by more than an order of magnitude. 

4 Outlook 

Feasible extensions in our framework include a calculation for the gauge cou- 
pling = 10, larger extensions in time or larger values of A. So far the 
LCE-programs were run on a SUN-workstation. Applied to the deconfincmcnt 
transition of QCD with dynamical fermions, the first approach would lead to a 
localization of the line of critical hopping parameters, whereas a derivation of 
an effective scalar model (in terms of quark condensates) along the lines of the 
second approach applies to the chiral transition of QCD. Finally we remark 
that DLCEs have a much wider range of applications than it was indicated 
here. They include spin glasses, partially annealed neural networks, or diluted 
Ising models and are independent of any variational approach. 
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